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We show that black phosphorus has room-temperature charge mobilities on the order of 10^ 
cm^V“^s“^, which are about one order of magnitude larger than silicon. We also demonstrate 
strong anisotropic transport in black phosphorus, where the mobilities along the armchair direction 
are about one order of magnitude larger than zigzag direction. A photocarrier lifetime as long as 
100 ps is also determined. These results illustrate that black phosphorus is a promising candidate 
for future electronic and optoelectronic applications. 


The rapid development of silicon-based integrated cir¬ 
cuit technology over the past several decades relied on 
shrinking the channel length of the transistors. This 
approach is facing a severe challenge posted by several 
fundamental limits. Hence, recently, significant efforts 
have been devoted to finding new materials for post¬ 
silicon electronic technology. One promising candidate, 
graphene, has shown room-temperature charge mobilities 
two orders of magnitude higher than siliconii^. How¬ 
ever, the lack of a band gap limits its application as 
a channel material in logic devices. Since 2012, tran¬ 
sition metal dichalcogenides, such as M 0 S 2 , have been 
extensively studied^. These monolayer semiconductors 
have sizable band gaps^ and possess several novel op¬ 
tical and spin/valley-related properties^. However, their 
charge mobilities are still lower than desired. 

Since 2014, two-dimensional black phosphorus (BP), 
also known as phosphorene, has emerged as a new 
promising material to provide both a sizable band gap 
and high charge mobilities. Few-layer and monolayer 
BP can be fabricated from bulk crystals by mechanical^ 
and liquid^ exfoliation. Although bulk BP has a rela¬ 
tively small band gap^ of 0.3 eV, several studies^dSi^— 
have shown that the band gap of BP thin films de¬ 
pends strongly on the thickness when approaching 
the monolayer limit - a feature similar to other two- 
dimensional materials: Charge transfer characteristics of 
field-effect transistors indicated a band gap of 1.0 eV in 
monolayersiS; while a 2.0-eV monolayer band gap was 
suggested by a scanning tunneling microscopy measure¬ 
ment on the surface of a bulk BP—. Photolumines¬ 
cence peaks from excitons were observed at 1.45 eV in 
monolayer BP—, and evolves from 1.3 to 0.9 eV when 
the number of atomic layers are changed from two to 
fivei^. Without knowing the exciton binding energies, 
these numbers give the lower limits of the band gap. Al¬ 
though more studies are still necessary to precisely de¬ 
termine the band gaps, the fact that few-layer BP pos¬ 
sesses band gaps large enough for logic applications has 
been established. On transport properties, theoretical 
studies indicated that monolayer and few-layer BP pos¬ 
sess high room-temperature mobilities'^ on the order of 
10 "^ cm^V“^s“^ and can outperform silicon and M 0 S 2 in 


ballistic devices^^^. Experimentally, field-effect transis¬ 
tors of multilayer BP have been fabricated^ii^^— . These 
transistors can have long operation lifetime once passi¬ 
vated by top AI 2 O 3 Iayers2ii22, and can operate at GHz 
regime^^. In optoelectronic applications, high photo- 
responsitivity^l and broadband response to both visible 
and infrared radiations^^ have been demonstrated in BP. 
Utilizing the ambipolar transport behavior, photovoltaic 
effect was demonstrated in few-layer BP in a local-gate 
p-n junction configuration^^. 

Here we report experimental observations of excep¬ 
tional transport performance of multilayer BP and its 
high in-plane anisotropy. By resolving spatiotemporal 
dynamics of photocarriers that are instantaneously in¬ 
jected by a tightly focused ultrashort laser pulse, we ob¬ 
tain room-temperature diffusion coefficients of 1,300 and 
80 cm^s“^ for photocarriers moving along the armchair 
and zigzag directions, respectively. These values corre¬ 
spond to ambipolar mobilities on the order of 50,000 and 
3,000 cm^V—ls“^, respectively. We also find that the 
photocarrier lifetime in multilayer BP is about 100 ps. 
These results provide fundamental parameters of BP, and 
illustrate its great promise for electronic and optoelec¬ 
tronic applications. 

Bulk BP crystal has orthorhombic symmetry, as il¬ 
lustrated in Fig. la, where single layers of phosphorus 
atoms are stacked by van der Waals force. In each layer, 
phosphorus atoms form a puckered honeycomb lattice, as 
shown in Fig. lb. Thin films of BP were mechanically 
exfoliated from bulk crystals purchased from 2D Semi¬ 
conductors. An adhesive tape is used to remove layers 
of BP from the crystal, which are then transferred to a 
Si 02 (90 nm)/Si substrate. The identified flake is imme¬ 
diately covered by a thin boron nitride layer, which was 
prepared with the same procedure, in order to prevent 
its degradation in air. The sample used in this study is 
shown in Fig. Ic. The thicknesses of the BP film and 
the boron nitride layer are determined to be 16 and 4 
nm, respectively, according to atomic force microscopy 
measurements. All the measurements were performed in 
ambient condition. No sample degradation was observed 
during the entire course of the study. 

In our transient absorption measurements, a 730-nm 
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FIG. 1. Lattice structure of black phosphorus and 
sample image, a, Lattice structure of multilayer black phos¬ 
phorus. b, Top view of a monolayer black phosphorus, with 
armchair and zigzag directions, c, Optical image of the sam¬ 
ple studied. The black phosphorus (BP) flake of 16-nm thick 
is partially covered by a 4-nm boron nitride (BN) flake, indi¬ 
cated by the dashed line. The horizontal (H) and vertical (V) 
directions in the laboratory frame are labeled. The angle 6 is 
between H and the long edge of the flake. 


and 100-fs laser pulse with a peak fluence of 10 /iJ cm“^ 
injects electron-hole pairs in the sample. The pulse is lin¬ 
early polarized along the long edge of the flake, as shown 
in Fig. Ic. The dynamics of the injected photocarriers 
are monitored by measuring the differential reflection of 
a 810-nm probe pulse that is polarized along the direc¬ 
tion perpendicular to the pump polarization. The dif¬ 
ferential reflection of the probe is defined as the normal¬ 
ized change in the probe reflection induced by the pump, 
AR/Rq = {R — Rn)/Ro- Here, R and Rq are the reflec¬ 
tion of the probe with the presence of the pump pulse or 
without it, respectively. As shown in Fig. 2a, the dif¬ 
ferential reflection signal reaches a peak near zero probe 
delay. After a transient process of about 20 ps, the signal 
decays exponentially (red line), with a time constant of 
100 ± 5 ps. We attribute this to the carrier lifetime in 
BP. By repeating the measurement with different values 
of the pump fluence (and hence different densities of the 
injected carriers), we confirm that the peak differential 
reflection signal is proportional to the pump fluence, as 
shown in the inset of Fig. 2a. 

The lattice structure of BP shown in Fig. la 
dictates that its in-plane properties are anisotropic, 
since the effective masses of the electrons and holes 
depend strongly on the crystalline directio n^^i^^~ — . 
Indeed, the anisotropic optical responses of BP 
have been predicte d^^i^° and experimentally confirmed 
by polarization-dependent Raman scatterin g^^d^ and 
absorption^. To study the anisotropic transient absorp¬ 
tion of BP, we measure the differential reflection signal 
with different pump-probe polarization configurations as 
we rotate the sample with respect to the polarization di- 
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FIG. 2. Time-resolved differential reflection signal 
from black phosphorus and its angular dependence. 

a. Differential reflection (AR/Rq) measured from the sample 
with a 730-nm pump and a 810-nm probe pulses. The red line 
indicates a single-exponential fit with a time constant of 100 ± 
5 ps. The inset shows linear relation between the peak signal 
and the pump fluence. b. The peak differential reflection sig¬ 
nal as a function of 6 , the angle between the long edge of the 
flake and horizontal direction in laboratory frame. The black 
squares, red circles, and blue triangles are measured with the 
pump-probe polarization configurations of horizontal-vertical 
(H-V), horizontal-horizontal (H-H), and vertical-vertical (V- 
V), respectively. 


rections. The squares, circles, and triangles in Fig. 2b 
show the peak signal of differential reflection as a function 
of the angle (9) between the long edge of the flake and the 
horizontal direction in laboratory frame (Fig. Ic), with 
the pump-probe polarizations of horizontal-vertical (H- 
V), H-H, and V-V, respectively. Pronounced anisotropy 
is clearly seen in all the configurations. 

The angular dependence of the differential reflection 
signal shown in Fig. 2b can be induced by both the pump 
absorption and the probe detection efficiency. First, due 
to the different effective masses along the armchair and 
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FIG. 3. Photocarrier diffusion along the armchair direction of black phosphorus, a, Differential reflection as a 
function of the probe delay and probe position obtained by scanning the probe spot with respect to the pump spot along the 
armchair direction for each probe delay, b, Examples of the spatial profiles of the differential reflection signal, obtained with 
probe delays of 2 (squares), 10 (circles) and 20 ps (triangles), respectively. The solid lines are Gaussian fits, c. The squared 
width as a function of the probe delay. The solid line is a linear fit that gives a diffusion coefficient of 1,300 ± 30 cm^s“^. 


zigzag directions, the absorption coefficients are differ¬ 
ent for pump pulses polarized along the two directions. 
Hence, as the sample is rotated with respect to the pump 
polarization, the pump injects different carrier densi¬ 
ties, which results in an angular dependent signal level. 
Second, the anisotropic bandstructure causes the probe 
pulses with different polarizations to sense the carriers 
with different efficiencies. Since the angular dependence 
in the H-V configuration is closer to H-H than V-V, we 
can conclude that the anisotropic differential reflection is 
mainly caused by the effect of pump absorption. When 
the pump polarization is along the long edge of the sam¬ 
ple, a maximal carrier density is injected. Hence, we 
assign the long edge to the armchair direction, which is 
expected to have the maximal absorptionii. At 0 = 0, 
the signal observed in the H-V configuration is about 
1.7 times larger than H-H, indicating that the V-probe 
is more efficient in detecting carriers than the H-probe. 
That is, the maximal probe efficiency is achieved when 
the probe polarization is along the zigzag direction. 

The solid lines in Fig. 2b are fits to the data by 
Acos^{9 + 0 q) + B. In the H-V configuration, we ob¬ 
tain a degree of anisotropy of {A + B)/B = 6.7. Given 
that the probe efficiency anisotropy is about 1.7, we de¬ 
duce a pump absorption anisotropy of about a factor of 
4. The pattern of H-H configuration is rotated with re¬ 
spect to H-V by about 8°, caused by a compromise of 
anisotropic pump and probe effects. The V-V pattern is 
rotated by 90° from H-H, which confirms that the angu¬ 


lar dependence observed is not from any artifacts of the 
experimental system. 

We study transport of photocarriers by performing 
spatially resolved differential reflection measurements. 
The tightly focused 730-nm pump injects photocarriers 
with a Gaussian spatial profile of about w = 2.4 /rm in 
full width at half maximum. Once injected, the photocar¬ 
riers diffuse in each atomic layer, driven by the in-plane 
density gradient. Hence, the spatial density profile ex¬ 
pands. This diffusion process is monitored by measuring 
differential reflection of the probe as a function of the 
probe delay and probe position. 

First, we study diffusion of photocarriers along the long 
edge direction of the flake, which is identified as the arm¬ 
chair direction according to Fig. 2b. In this measure¬ 
ment, both pump and probe polarizations are along the 
armchair direction, and the probe spot is scanned along 
the same direction. Figure 3a shows the differential re¬ 
flection signal as a function of the probe delay and probe 
spot position with respect to the pump spot. At each 
probe delay, the profile has a Gaussian shape. Figure 3b 
shows three examples of these profiles at probe delays of 
2 (squares), 10 (circles), and 20 ps (triangles), respec¬ 
tively, along with the corresponding Gaussian fits (solid 
lines). From each fit, we obtain the width (full width 
at half maximum) of the profile. The squared width is 
plotted in Fig. 3c as a function of the probe delay. 

According to the diffusion model^, the profile of car¬ 
rier density remains its Gaussian shape during the diffu- 
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FIG. 4. Photocarrier diffusion along the zigzag direction of black phosphorus, a, Differential reflection as a function 
of the probe delay and probe position obtained by scanning the probe spot with respect to the pump spot along the zigzag 
direction for each probe delay, b, Examples of the spatial profiles of the differential reflection signal, obtained with probe 
delays of 2 (squares), 10 (circles) and 20 ps (triangles), respectively. The solid lines are Gaussian fits, c. The squared width as 
a function of the probe delay. The solid line is a linear fit that gives a diffusion coefficient of 80 ± 5 cm^s“^. 


sion, as confirmed in Fig. 3b. The width of the profile 
evolves as = w^{to) + 161n(2)£)(t — to), where to is 
an arbitrarily chosen initial time, and D is the diffusion 
coefficient. Such a linear expansion of the area covered 
by the photocarriers is confirmed by the linear fit to the 
data shown in Fig. 3c (red line). From the slope, we 
obtain a diffusion coefficient of 1,300 ± 30 cm^s“^. 

Next, we rotate the sample by 90° and repeat the mea¬ 
surement under otherwise the same conditions. Now the 
polarizations and the scan directions are all along the 
zigzag direction. The results are presented in Fig. 4, 
in the same fashion as Fig. 3. Clearly, the diffusion is 
much slower along the zigzag direction. From a linear fit, 
shown as the solid line in Fig 4c, we obtain a diffusion 
coefficient of 80 ± 5 cm^s“^, which is about a factor of 
16 smaller than along the armchair direction. 

The spatially resolved measurements also help clarify 
the fast transient process observed in Fig. 2a, where 
the signal dropped by about a factor of 4 in the first 
20 ps. As a rough estimation from Figs. 3c and 4c, 
the area of the profile increases by more than a factor 
of 2. Hence, the peak of the profile should decrease by 
the same factor due to the diffusion along. On the other 
hand, at later probe delays, the diffusion makes smaller 
contribution to the decay of the signal since the profile is 
already significantly broadened. This is confirmed by the 
single-exponential decay as indicated by the red line in 
Fig. 2a and validates the attribution of the 100-ps time 
constant to the carrier lifetime. 

Our noninvasive and all-optical measurement of spa- 


tiotemporal dynamics of photocarriers in BP reveals its 
exceptional and highly anisotropic transport properties, 
which show the great potential of this material for elec¬ 
tronic and optoelectronic applications. The measured 
diffusion coefficient can be used to deduce the mobility 
by using the Einstein’s relation, /i/e = D/ksT, where /i, 
e, kB, and T are the mobility, elementary charge, Boltz¬ 
mann constant, and temperature. We obtain mobilities 
of 50,000 ± 2,000 and 3,000 ± 200 cm^V-is-^ along 
armchair and zigzag directions, respectively. 

In the diffusion process we studied, the photocarriers 
move as electron-hole pairs due to the Coulomb attrac¬ 
tion. Hence, the diffusion coefficient obtained are am- 
bipolar diffusion coefficient, which is related to the dif¬ 
fusion coefficients of electrons and holes {Df. and Dh) by 
Da = 2DeDh/{De + Dh) for photocarriers^. In gen¬ 
eral, the ambipolar diffusion coefficient is controlled by 
the slower carriers. Hence, the ambipolar diffusion coeffi¬ 
cient can be regarded as the lower limit of D^ and Dh ■ In 
BP, the effective masses of electrons and holes are within 
a factor of two^. Hence, the ambipolar diffusion coeffi¬ 
cient and ambipolar mobility provide satisfactory order- 
of-magnitude estimation of the diffusion coefficients and 
mobilities of both electrons and holes. 

Recently, high room-temperature charge mobilities on 
the order of 10^ cm^V“^s“^ in BP has been predicted 
theoretically^^. The previously reported experimental 
values have yet to reach that level. For example, transis¬ 
tors of 10-nm BP flakes achieved mobilities in the range 
of 100 - 1000 cm^V“^s“^ at room temperature^ d^d^ , 
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and higher values at low temperature a^^d^ . A trilayer 
BP transistor demonstrated electron and hole mobilities 
of 38 and 172 cm^V“^s“^, respectively^^. The electri¬ 
cal measurements may not reveal the intrinsic transport 
properties of the material since the transport can be lim¬ 
ited by grain boundaries and by the quality of the elec¬ 
trodes. Indeed, significant influence of the metal contact 
has been suggested^^— . On the other hand, the all- 
optical approach we used in this study does not involve 
device fabrication and grain boundaries play a negligible 
role in the diffusion process. We find that our results are 
reasonably consistent with the theoretical predictionsi^. 

The anisotropic transport property can be mainly 
attributed to the anisotropic effective masses of elec¬ 
trons and holes. The effective masses of electrons with 
wavevectors along the armchair and zigzag directions are 
0.0825mo and 1.027mo, respectively, where mo stands 
for the electron rest mass in vacuum. Hole effective 


masses are 0.0761mo and 0.648mo along the two direc¬ 
tions, respectively^. The differences in effective masses 
suggest that the mobilities along the armchair direction 
are about one order of magnitude larger than zigzag 
direction. Such a ratio is also predicted in monolayer 
and few-layer Based on our result, the mobil¬ 

ity along the armchair direction is about 16 times larger 
than the zigzag direction, which is reasonably consistent 
with these predictions. 
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